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We analyze a lattice model closely related to the one-dimensional inelastic gas with periodic 
boundary condition. The one-dimensional inelastic gas tends to form high density clusters of par- 
ticles with almost the same velocity, separated by regions of low density; plotted as a function of 
particle indices, the velocities of the gas particles exhibit sharp gradients, which we call shocks. 
Shocks and clusters are seen to form in the lattice model too, although no true positions of the 
particles are taken into account. The locations of the shocks in terms of the particle index show re- 
markable independence on the coefficient of restitution and the sequence of collisions used to update 
the system, but they do depend on the initial configuration of the particle velocities. We explain 
the microscopic origin of the shocks. We show that dynamics of the velocity profile inside a cluster 
satisfies a simple continuum equation, thereby allowing us to study cluster-cluster interactions at 
late times. 

PACS numbers: 47.70.Nd, 45.70.Mg, 05.40.-a, 81.05.Rm 



I. INTRODUCTION 

Dynamics of granular fluids has captured a lot of at- 
tention from theoretical physicists for the last few years. 
In a theorist's model, the constituent particles of a gran- 
ular fluid, usually considered to be hard spheres of finite 
radii, irretrievably dissipate kinetic energy via inelastic 
binary collisions and interparticle frictional forces. As a 
result, unlike the microscopic models for the classical ki- 
netic theory of g granular fluid that is not driven 
by external forces "cools freely" . 

Even in the absence of frictional forces, in a stark con- 
trast to hard sphere fluids with elastic inter-particle bi- 
nary collisions, such a simplified model of granular flu- 
ids exhibits complex behaviour at macroscopic scales. In 
its simplest form, a freely cooling initially homogeneous 
and isotropic (both in the particles' position and velocity 
space) inelastic gas spontaneously forms nontrivial struc- 
tures in the macroscopic velocity as well as in the macro- 
scopic density field of the gas. A large number of studies, 
mostly from the point of view of inelastic hydrodynam- 
ics, have been carried out to understand the onset for the 
formation of these structures in two and three dimensions 
(Refs. QESSHEEIto cite a few). The (qualitative 
and quantitative) picture that has emerged from these 
studies is that in two and three dimensions, the system 
of inelastic hard spheres suffers from inherent long wave- 
length linear instabilities. When the system size allows 
such long wavelengths to be present, these instabilities 
start to generate inhomogeneities in the macroscopic ve- 
locity and the density field of the gas. In the subsequent 
evolution, these inhomogeneities interact nonlinearly to 
give rise to macroscopic structures and the entire sys- 
tem evolves into a collection of densely populated clus- 
ters that are separated by regions containing particles at 
very low density [lj,|8|. At late times, the clusters collide 
in a very complex manner and merge — a phenomenon 
that is known as coarsening in the literature 0, ^| . 



The late time evolution of a freely cooling inelastic gas 
is thus qualitatively completely different from the linear 
instability mechanisms at early times. However, although 
a very large number of studies have been devoted to ki- 
netic theory of freely cooling inelastic gases in two and 
three dimensions, due to the difficulties associated with 
the nonlinearities in the behaviour of individual clusters 
and cluster-cluster collisions, a proper theoretical under- 
standing of the long time dynamics of freely cooling in- 
elastic gases has remained elusive. The existing results 
have only been numerical 0, U [H GH 13 El 

At the other extreme, fully analytical solutions have 
been found for completely inelastic (or "stic ky" ) granu- 
lar gas (of point particles) in one dimension [lj], and it 
has been shown that the sticky gas in one dimension is 
described by the Burgers equation in the inviscid limit 
In addition, a recent experiment has also observed 
clustering in a one-dimensional granular gas ^3 . Due to 
dimensional reasons, the dynamics of a granular gas in 
one dimension is qualitatively different from those in two 
or three dimensions (e.g., vortices cannot form, a strict 
ordering of particles from left to right is maintained at all 
times), but structures are still seen to form in the veloc- 
ity and as well as in the density field of the gas 0, . 
Ben-Nairn et al 17J have studied the formation of these 
structures, and have conjectured that at the long times, 
the behaviour of a one-dimensional inelastic gas should 
be the same as that of the sticky gas. 

Our purpose in this paper is to unravel some of these 
long-time phenomena at a more microscopic level in a 
simple one-dimensional lattice model that has been in- 
troduced in Ref. [Til Il9|. In this model, one considers 
a system of N particles on the integral lattice positions 
(denoted by k) of a ring of size N with initial velocities 
chosen randomly from a uniform distribution in [—1, 1]. 
The ordering of the particles is maintained at all times; 
however, since the particles do not move in this lattice 
model, the velocities of the particles are not the time 
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derivative of their positions. Instead, Vk, the velocity of 
the fcth particle, is simply a scalar quantity associated 
with the fcth particle, and Vk changes only when the fcth 
particle participates in a collision with one of its neigh- 
bours according to the following collision rule. At a col- 
lision between kth and (k + l)th particles at any time, 

the post-collisional velocities ^ k + (k+i) are re l a ted to their 



pre-collisional velocities v 
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where < r < 1 is the coefficient of restitution. Time is 
measured by the average number of collisions per particle 
in this lattice model. 




FIG. 1: Comparison of the shock locations on a ring of size 
N as a function of particle positions k, k — 1, 2, . . . , N: (a) 
a snapshot of one- dimensional inelastic gas |17| for r — 0.3, 
(b) a similar snapshot of the random lattice model for r = 
0,(c)a snapshot of the systematic lattice model. All systems 
had identical initial positions and velocities of the particles. 
Visual inspection shows that the locations of the (relative) 
large shocks are aligned almost perfectly. The scale in the 
y-direction is arbitrary. 



lision frequency is chosen to be independent of the in- 
coming velocities of colliding particles. For that model, 
it has been analytically shown that correlations develop 
with a diffusively growing correlation length, which con- 
sequently affects the temperature decay rate. In addition, 
the inelastic lattice model in one-dimension with the ki- 
netic constraint, which we study here, was also analyzed 
in Refs. |18l llP| in terms of velocity distribution and 
structure factors |2l|. It is however important to real- 
ize the difference between these existing results and the 
ones reported in this paper: the existing results mainly 
concern global quantities, while in this paper, our main 
thrust is to study the behaviour of the microscopic inho- 
mogeneities arising from the kinetic constraint. 

The key feature of both variants of the lattice model 
with kinetic constraint studied here is that an initial con- 
figuration of random velocities of the particles soon de- 
velops distinct spatial structures, eventually leading to 
large positive sharp jumps in the particle velocities in- 
terspaced with relatively smooth variations. We refer to 
the large positive sharp jumps in the velocity field as 
shocks and the region between two consecutive shocks, 
where the velocity variations are smooth, as clusters. 
These structures have already been observed in , but 
have not been fully analyzed. Formation of shocks and 
the subsequent dynamics of the clusters make the lattice 
model interesting on its own, but its relevance is realized 
only when the locations of its shocks in w(fc)-profile are 
compared to those of the one-dimensional inelastic gas 
[3 113 (see Fig. HJ). By contrast, there are no shocks 
in the one dimensional random lattice model without ki- 
netic constraints [2p| . 

This paper is organized in the following manner: in 
Sec. |nl we discuss the generic features of the lat- 
tice model and explore its connections with the one- 
dimensional inelastic gas [Til Il7j . In Sec. IIIII we an- 
alyze the formation of shocks at early times and cluster 
dynamics at late times. We finally end the paper with a 
short discussion in Sec. II VI 



II. GENERIC FEATURES OF THE LATTICE 
MODEL 



In this paper, we will consider two variants of this lat- 
tice model, namely the "random lattice model" and the 
"systematic lattice model" . Both obey the collision rule 
JTJ , but they differ in the way a colliding pair of particles 
is chosen to update the system. For the random lattice 
model, at any update, the colliding pair of particles are 
chosen randomly from all particles that momentarily sat- 
isfy the kinetic constraint Vk+i — Vk < 0. On the other 
hand, for the systematic lattice model, the colliding pair 
of particles are the ones that has the momentary global 
minimum value of Vk+i — Vk- 

The random lattice model without kinetic constraint, 
has already been studied as part of a larger class of mod- 
els, the inelastic Maxwell models [2(J, in which the col- 



A. Phenomenology of shock development and the 
subsequent dynamics for the random lattice model 

To start with, in Fig. [3 we show a time sequence of 
the velocity profile v(k) for the random lattice model of 
10000 particles with r = 0.7. The initial configuration 
[Fig. HJa)] is created by choosing Vk randomly from a 
uniform distribution in [—1,1]. As can be seen in Fig. 
|2{b), velocity correlations set in very rapidly (within 10 
collisions per particle). After 10 4 collisions per parti- 
cle, shocks and clusters can be clearly identified [Fig. 
Hfc)]. According to the kinetic constraint, neighbouring 
particles cannot collide across a shock, and as a result, 
each cluster evolves independently of the others, until 
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two neighbouring clusters collide and coalesce to form a 
new bigger cluster [this mechanism is illustrated by the 
evolution of velocity profile of Fig. Etc) to that of Eld); 
Fig. Eld) corresponds to a state after 10 6 collisions per 
particle] . 

From the above phenomenological description, one can 
identify two separate regimes: the initial homogeneous 
regime, and the subsequent clustered regime at late 
times. To make the evolution of the system from one 
regime towards the other quantitatively more precise, we 
notice that it corresponds to a symmetry breaking of the 
probability distribution P(Avk,t) of relative velocities 
Awj, = — Wfe. Initially, for each fc, the distribution of 
Vk is uniform in [— 1, 1] (i.e., a "box" distribution). Nat- 
urally, P(Avk, t = 0) is a piecewise linear distribution on 
[—2,2] (a "triangular" distribution), symmetric around 
Avk — 0. Under the effect of the dynamics, this sym- 
metry is not preserved. This is easily understood from 
the fact that only the neighbouring particle pairs with 
Avk < can collide, and after a collision between the kth 
and (fc + l)th particle, Avk-i, Avk and Avk+i change in 
a way that, in general, does not preserve the symmetry of 
P(Avk,t). As this process continues, from phenomeno- 
logical considerations, we know that a few shocks remain 
at late times. Now, since shocks correspond to large val- 
ues of Avk, this means that in the clustered regime, the 
probability to have a large positive value of Avk is ex- 
pected to be greater than that of a correspondingly large 
negative value. 

To detect the locations of the shocks by following the 
evolution of P(Avk,t) as a function of time, one can 
identify the location k of a shock by requiring Avk > 
C Imim, Avj\, where C > is a constant. The choice of 
the numerical value of C is arbitrary, and this arbitrari- 
ness can be used to tune the minimum value that Avk 
must have in order to qualify for a shock. Furthermore, 
if one defines the instant t p when the system passes from 
one regime to the other by the minimal time where at 
least one shock becomes visible, then the above require- 
ment can also be used to characterize t p . Of course, the 
precise value of t p depends on the chosen value of C. 

In the clustered regime, due to inelastic collisions, the 
amplitude of the velocity profile within each cluster de- 
creases, and the velocities of the particles approach that 
of the center of mass of the cluster itself. This process 
continues for a while until two neighbouring clusters, 
with center of mass velocities Vj and V^+i collide and 
coalesce and the shock separating them disappears [at 
time t, if we number the clusters j = 1,2, . . . ,ji(t) from 
the left to the right, then such a cluster-cluster collision 
takes place only if Vj > Vj+i]. Precisely this mechanism 
is responsible for making the velocity profile of Fig. Et c ) 
evolve to that of Fig. Eld). 

At very late times, due to periodic boundary condi- 
tion, only a single cluster survives. Its amplitude also de- 
creases in time, and the velocities of all particles asymp- 
totically converge to the center of mass velocity of the 
whole system. 




FIG. 2: An example time sequence of shock developments for 
a system of 10000 particles in the random lattice model: (a) 
Initial profile (randomly chosen Vk from a uniform distribution 
in [—1, 1]), (b) after 10 collisions per particle, (c) after 10 4 
collisions per particle, (d) after 10 6 collisions per particle. 
The largest \vk ± — vi C2 \ Vfci,fc2 have been scaled to unity in 
each graph. 



B. Dependence on the coefficient of restitution for 
the random lattice model 

One of the first questions that one can ask for the lat- 
tice model is what is the effect of inelasticity on the clus- 
tering properties of the system. Clearly, for r = 1, the 
collisions are elastic, and upon collision, the velocities of 
two neighbouring particles are simply interchanged. Such 
a dynamics cannot lead to clustering or any structure for- 
mation. 

On the other hand, our simulations show that the lat- 
tice model exhibits clustering for any value of r smaller 
then 1, and it is only t p defined in Sec. Ill Al that be- 
comes an increasing function of r. The most remarkable 
observation, however, is that for a given initial velocity 
configuration of the particles, not only the locations of 
the shocks persisting at late times are the same for all r, 
but also the macroscopic velocity profiles are nearly in- 
dependent of r. If we consider the macroscopic velocity 
profile of the lattice model with a coefficient of restitution 
n at a time t\ well into the clustered regime, then for 
any r2 > r\, (originating from the same initial velocity 
configuration of the particles) there exists a time ti < t\ , 
for which the two macroscopic velocity profiles are almost 
the same (see Fig. EJl. Effectively, this implies that for 
all r, the dynamical behaviour of the system is identical 
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to that for r = at late times. We also observe that at 
intermediate times, the maximum number of shocks (and 
correspondingly, clusters) in the system are observed for 
r = (see Sec. IIII Al in this regard). 




FIG. 3: Velocity profiles for the random lattice model of 5000 
particles: (a) for r — 0.7 after 10 4 collisions per particle (b) for 
r — 0.95 and 10 3 collisions per particle. The initial velocity 
configurations of the particles were the same for both (a) and 
(b). Once again, the largest \vk 1 —Vk 2 1 Vfci, fe has been scaled 
to unity for both profiles. 



colliding pair is automatically determined by the instan- 
taneous minimal value of (xk+i — %k)/( v k ~ Vk+i)j for 
Vk — i'k+i > 0. In this sense, one can view the one- 
dimensional inelastic gas as a variant of the lattice model 
too, where the colliding particles are chosen in a very 
complicated manner. Section Hi CI then indicates that for 
a given initial configuration of the particles' velocities 
as a function of the particle index, the locations of the 
shocks in the one-dimensional inelastic gas and those in 
the lattice model are the same. That this is indeed the 
case is numerically verified in Fig. ^ 

In our simulations of the one-dimensional inelastic gas, 
we followed the procedure outlined in Ref. ^7|> namely 
that to avoid the inelastic collapse, the collisions were 
taken to be elastic when the relative velocities between 
the neighbouring particles are smaller than a given cutoff. 
Our results confirm the existence of "Burgers" shocks in 
the v(x) profile ^tJ- More importantly, we note that 
for the one-dimensional inelastic gas, there is an inverse 
correspondence of shocks and clusters between this v(x) 
and v(k) profile — each shock in the v{x) profile [i.e., a 
large particle density n(x)] corresponds to a cluster in the 
v(k) profile and the shocks in the v(k) profile correspond 
to regions with finite gradients in v(x). 



C. Dependence on the order of collisions and on 
the initial velocity configuration of the particles 

The next important question one can ask is if not r, 
exactly what influences the positions of the shocks? A 
priori, one expects that the locations of the shocks and 
the long time dynamics of the system depend on the ini- 
tial velocity configuration of the particles, as well as on 
the sequence in which the collisions are performed. 

Surprisingly however, the locations of the shocks de- 
pend only on the initial configuration of the particles' 
velocities, e.g, two different random collision sequences 
result in the same locations of the shocks. On the other 
hand, the variant of the lattice model (random or sys- 
tematic) determines the macroscopic shape of the veloc- 
ity profile inside the clusters. To demonstrate this phe- 
nomenon, starting from the same initial velocities of the 
particles, we simulated both the random and the sys- 
tematic lattice models and compared the locations of the 
shocks (see Fig. P). 



D. Relation between the lattice model and 
one-dimensional inelastic gas 

The empirical observations in the preceding sections 
suggest that so long as one is only interested in the ve- 
locity profile v(k), a close relation exists between the 
one-dimensional inelastic gas and the lattice model from 
the following consideration: the dynamics of a one- 
dimensional inelastic gas for a given initial configuration 
is completely deterministic — after any collision, the next 



III. ANALYSIS 

A. Development of shocks 

The lesson that we learnt from Sec. [H] is that the 
locations of the shocks in this lattice model are essentially 
determined from the initial configuration of the particles' 
velocities. This immediately gives rise to the following 
question: how can one predict the locations of the shocks 
seen at early times from the initial velocity configuration 
of the particles? 



v(k) 




FIG. 4: A schematic diagram of the initial velocity profile for 
a few particles to explain the development of shocks. 



As it turns out, most of the shocks observed in the 
clustered regime already exist at t = 0. These are the 



locations marked by large positive jumps in the v(k) pro- 
file between two neighbouring lattice sites at t = 0, sur- 
rounded by relatively small velocity variations. To illus- 
trate this point, we present a schematic diagram for such 
a velocity profile for a few particle in Fig. In the 
subsequent dynamics at short times, the particles on the 
right of B and on the left of A collide with each other and 
"thermalize" (i.e., the velocity variations reduce due to 
inelasticity), but this thermalization may not necessarily 
convert a relatively large vb — va > to vg — < 
quickly enough to make the A and B collide. In other 
words, such a large positive jump remains preserved and 
may eventually give rise to a shock. 

To check if such a scenario is correct or not, one needs 
to identify the locations of such large positive jumps in 
the initial v(k) profile, and contrast them with the ob- 
served locations of the shocks at early times. A conve- 
nient way to single out such jumps [say, in Fig. Ufa)] 
is to define a (2n + l)-coarse grained velocity profile 

k+n 

v(k) = Vi/(2n + 1) of the particles, which, roughly 

i—k — n 

speaking, can single out one such large jump over the sur- 
rounding small velocity variations in a window of (2n+ 1) 
lattice sites in the v(k) profile. Having fine-tuned n, a 
surprisingly close match between the observed shock lo- 
cations at early times and the initial configuration of the 
particles' velocities can be found (see Fig. EJ). 

Notice, however, in Fig. 0that not all the large jumps 
in the initial configuration of the particles' velocities have 
turned out to be shocks at a later time. Indeed, in gen- 
eral, whether such an initial large positive jump at a 
particular location develops into a shock or not really 
depends on the magnitude of the jump itself in relation 
to particle velocities in its immediate vicinity, and as well 
as on the coefficient of restitution r. In this regard, the 
dependence on the coefficient of restitution is not difficult 
to understand, as r quantifies the "transport of velocity" 
between two colliding particles. Intuitively speaking, an 
inelastic collision between two neighbouring particles can 
be viewed as a combination of dissipation (thermaliza- 
tion) and transport of velocities. For r — 1 only trans- 
port of velocities can take place. At the other limit, for 
r = 0, there is no transport of velocities but only dissi- 
pation. As a result, for r — 0, small velocity variations 
around any large jump thermalize immediately, leading 
to an early appearance of the clustered regime, and a 
large number of initial jumps end up becoming shocks. 
With increasing r, some of the initial large positive jumps 
are eliminated by transport of velocities, and clustering 
regime, with a smaller number of selected shocks (and 
clusters), appears later. 

The appearance of the shocks at the locations where 
relatively large positive jumps exist is therefore caused 
by the thermalization of the smaller velocity variations 
around these jumps — due to thermalization, the smaller 
velocity variations become even smaller, but the large 
positive jumps remain preserved under the dynamics. In 
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FIG. 5: Correspondence between the coarse grained velocity 
profiles v(k) (n — 7) for r — for the random lattice model of 
1000 particles. Gray curve corresponds to t = and the black 
curve represents early shocks. The scale in the y-direction is 
arbitrary. 



relation to their surroundings, the magnitude of these 
large jumps thus starts to grow in time. Such a scenario 
applies independently of the details of the model, and 
gives a qualitative explanation for the observation that 
starting from a given initial configuration, shocks appear 
at the same locations in the v(k) profile for both variants 
of the lattice model and as well as for the one-dimensional 
inelastic gas jQ, ^J. Moreover, it also clearly demon- 
strates that the shock formation "instability" in the v(k) 
profile in all these models is a "relative instability" and 
not an "absolute" one. Such behaviour has also been 
found in the (long-wavelength) linear instability of the 
macroscopic flow field for two- and three-dimensional in- 
elastic gases (see Refs. 0,0,0 for example). 



B. Cluster dynamics 

1. Dynamics within a single cluster 

Once the shocks develop, our simulations show that for 
a given variant of the lattice model, the velocity profile 
within each cluster has the same characteristic macro- 
scopic shape. For the deterministic lattice model, this 
shape is simply linear with a negative slope. Although 
we cannot provide an analytical derivation of this shape, 
if we interpret the dynamics along the lines that at each 
updating step of a cluster, two neighbouring particles 
collide at the location where the Vk+\ — Vk is minimum 
(which is a mechanism that tries to align all Vk+i — «fe 
values to its maximum value within each cluster), then 
the linear profile appears to be intuitively reasonable. 

For the random lattice model, however, the slope of the 
macroscopic velocity profile is smoothly varying, as can 
be seen in, e.g., Fig. |2fd) [occasionally, flat profiles can 
be observed too, but for the time being, we leave them 
for Sec. IIII B 2) . It turns out that the functional form of 
the smooth macroscopic velocity profile can be obtained 
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through a mean-field approach by averaging over an en- 
semble of random collision sequence realizations within a 
cluster, as we describe below. 

The idea behind this mean-field approach is the fol- 
lowing: we consider a cluster of M particles, and denote 
the velocity of the fcth particle at time t, averaged over 
random collision sequences, by (vk(t)). From the very 
definition of a cluster in the introduction, we assume 
that {vk+i{t)) — (vk(t)) < at all times. Thereafter, 
as we observe from direct simulation measurements that 
the probability of collisions between any two neighbour- 
ing particle pairs are equally likely, in this mean field 
approach, (vk(t)) is easily seen to satisfy the equation 

(« fc (t + l)> = (l-2p)(« fc (t)) 

+ p[(« fc (t)> - e{(v k (t)) - (« fc -i(t))}] 
+ p[(v k (t)) + s{(v k+1 (t)) - (v k (t))}]-m 

Here, p is an effective probability of a collision between 
any two neighbouring particles in this mean-field theory 
and e = (1 + r)/2. The three terms on the r.h.s. of 
Eq. (JjjJ respectively originate from the events when the 
fcth particle is not involved in a collision, when there is 
a collision between the fcth and the (k — l)th particle 
and when there is a collision between the fcth and the 
(k + l)th particle. As Eq. (j2J) very simply reduces to 

(v k (t + 1)) = (1 - 2pe) (v k (t)) + pe <w fc _i(t)> 

+ pe(v k+1 (t)), (3) 

the interesting point to note is that Eq. is the discrete 
(both in space and time) form of the diffusion equation 

8{v(k,t)) 



at 



= D(r) 



d 2 {v(k,t)) 
dk 2 



(4) 



where mean- field theory predicts D (r) = pe. One has to 
keep in mind however that Eq. ^j) holds only for mono- 
tonically decreasing (v(k,t)) as a function ofk. A similar 
equation has also been found in the one-dimensional lat- 
tice model without the kinetic constraint [20| . 

To solve Eq. (0}, the macroscopic velocity profile of a 
cluster at time to can be generally expanded in a Fourier 
series. Since particle velocities are not transported across 

the boundaries of a cluster, d^fcft)) mus ^ van i s h at the 

ok 

boundaries of a cluster, the Fourier series contains only 
terms oc cos[tt jk]. The amplitude of the jth such term 
decreases as e _7r J D(r){t-t ) ^ SQ f or i ar g e va i ue s of 
(t — to), only the slowest decay mode j = 1 survives. 
Hence at long times we expect (v{k,t)) to be given by 



irk 
~M 



(v(k,t)) = A(t) cos 

= A(to)e-^ D(r ^-^/ M ' A 



nk 
~M 



(5) 



The mean-field result (JSJ) allows us to compare the half- 
cosine form of (v(k, t)} with the macroscopic velocity pro- 
file of the particles within a cluster, as observed in the 
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FIG. 6: (a) Comparison of the half-cosine solution 10 [solid 
line] and simulation data[open circles], the two are completely 
indistinguishable from each other, (b) Numerically obtained 
dependence of D(r) on r for the random lattice model. The 
simulations were carried out for M = 1000. The scale along 
the y axis for Fig. (a) is arbitrary. D(r) exhibits a singularity 
at r = 1, which is in clear disagreement with the mean-field 
theory prediction. 



simulation for the random lattice model. This compari- 
son is shown in Fig. Ela) — the fact that we cannot dis- 
tinguish the simulation data from the half-cosine shape 
of (v(k, t)) as predicted by Eq. JSJ) is an indication of how 
well the mean-field approach works to describe the aver- 
age shape of an isolated cluster. In addition, by following 
the dynamics from the simulation, we can also obtain an 
empirical functional form of D(r) [shown in Fig. EJb)]. 

There are two features of the D(r) vs. r curve that 
require further elaboration. The first one of them is that 
D(r) is an increasing function of r, indicating that the 
amplitude of the cluster A(t) decreases faster as r in- 
creases. This observation seems to contradict the fact 
that the dissipation decreases with increasing r. The 
point to notice however is that the mean-field approach 
really describes the time evolution within a cluster in 

which ^ k( U. < o. From that point of view, the ampli- 
ok 

tude of the cluster decreases not only from the dissipa- 
tion, but also from transport of velocities. This is easily 
seen from the fact that after a collision, the ordering of 
the magnitudes of the velocities of the colliding particles 
are simply exchanged, i.e., the particle on the right has 
a higher velocity, which results in a local positive slope 
in the v(k) profile. For increasing r, transport of veloc- 
ities become more effective, and thus the negative slope 

r d(v k (t)) 

of — — within a cluster decays to zero faster. Only 

if r = 1, the velocities of the particles within an isolated 
cluster order themselves in an increasing order of mag- 
nitude. However, D(r) has a singularity at r = 1 and 
the r = 1 case cannot be treated within the scope of 
our mean-field theory. Secondly, from the way we pre- 
sented our mean-field approach, it may seem that the 
r-dependence of D(r) can be obtained through the de- 
pendence of p and e on r. This is actually not true. 
The mean-field approach J2J neglects fluctuations in the 
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particles' velocities around (v(k, t)), but in an actual sim- 
ulation, these fluctuations are very important as they de- 
cide the sign of Av/.(t) — Vk+i(t) — Vk(t) for the particles 
within the cluster and thereby control which collisions 
are possible and which are not. 



2. Interacting clusters for the random lattice model 



With the background of Sec. IIII B II it is now clear that 
as time progresses, the amplitude of each cluster present 
in the random lattice model at late times effectively de- 
cays exponentially as ~ exp[— 7r 2 D(r)(t — to)/M 2 ]. Such 
a decay brings the velocities of the particles within a clus- 
ter closer and closer to the centre-of-mass velocity of the 
cluster itself. Let us assume that at time to 3> 1, there 
are j\ number of clusters present in the whole system of 
the random lattice model. If we number these clusters by 
j such that j = 1,2,..., ji and denote the number of par- 
ticles within the clusters, the amplitudes of the clusters 
and their velocities of their centre of masses respectively 
by Mj , Aj and Vj , then with increasing time, we observe 
the following dynamics: (i) when Vj-i < Vj < Vj+i, 
the shocks on the left and on the right of the jth clus- 
ter cannot respectively be smaller than Vj — Vj-i > 
and Vj+i — Vj > 0. In that case, the amplitude Aj sim- 
ply decreases to 0, forming a "flat" velocity profile (see 
e.g., Fig. |3J). (ii) On the other hand, when Vj > Vj + i , the 
magnitude of the shock between the jth and the (j + l)th 
clusters decreases to zero, and the two clusters coalesce 
together to form a new bigger cluster in a finite time. For 
a given configuration of clusters at time t, the mean-field 
theory of Sec. IIII B II provides us with a way to mea- 
sure the time t' = t + At m i n of the first coalescence of 
two clusters in the system. In fact, At m i n = min-,- Atj, 
where for each pair of neighbouring clusters (j, J + with 
Vj — Vj + i > 0, Atj is obtained by solving the equation 



A j+1 (t ) cxp[-n 2 D(r)At 3 /Mf +1 ] 

+ Aj(t ) exp[-TT 2 D(r) Atj/M?] = Vj - V J+1 . (6) 

Due to the fact that Eq. (jHJ) is transcendental in na- 
ture, a closed form analytical solution for Atj is impossi- 
ble to obtain, let alone the value of At m i n . Nevertheless, 
Eq. JBJ) provides us with a glimpse of how complicated it 
is to theoretically study the cluster-cluster collisions and 
coalescence in the random lattice model, and in general, 
in inelastic gases. If two clusters I and I + 1 are the first 
ones to coalesce (at time to + Ati) in the random lattice 
model, then a new cluster with Mi + M;+i particles and 
center of mass velocity (MiV + Mi +1 Vi+i) / {Mi + A//+i) 
is formed. At time to + Ati, the shape of the new cluster 
is different from a half-cosine, but as the mean-field the- 
ory of Sec. IIII B II suggests, very soon the shape of the 
new cluster converges to a half-cosine, unless the newly 
formed cluster collides with another one in the meantime. 



A "hydrodynamic" description of the random lattice 
model 



So far, we have analyzed the system-wide properties 
of shocks and clusters. We have also seen that there ex- 
ists an effective dynamics in terms of diffusion equations 
(0J and (0 within a single cluster. Based on this col- 
lected wisdom on the random lattice model so far, one 
can naturally ask if it is possible to express the system- 
wide properties of shocks and clusters in terms of an ef- 
fective ("hydrodynamic") equation. 

It turns out that indeed such an equation can be con- 
structed for the random lattice model: 



Or 

at 



D(r) [6(-V H 



e(-V - u) V~v] , (7) 



where denotes a unit step function, (\7 + v)k = v k+i — Vk 
and (V~t>)fc = Vk—Vk-i are the discrete gradients operat- 

d 

ing to the right and left respectively and — the discrete 

time gradient. It is important to stress that in a nu- 
merical implementation of Eq. Q , all pairs of particles 
with negative relative velocities collide simultaneously in 
a unit timestep. The main difference between the numer- 
ical implementation of Eq. Q and the random lattice 
model lies in the order of the collisions, but in the light 
of the previous results, we expect to find the shocks at 
the same locations. Moreover, inside a given cluster, 1(7)1 
reduces to (@J , so that the cluster velocity profiles should 
be the same as in the random lattice model. 




FIG. 7: Comparison of the macroscopic velocity profile (solid 
curve) predicted by Eq. J7J with the corresponding actual 
simulation data (open circles) of the random lattice model 
for r = 0. The simulations were carried out for 1000 particles 
and the largest \vk x — Vk 2 1 Vfci, fea have been scaled to unity 
both for the actual simulation data and the solid curve. 

Indeed, we find that the macroscopic velocity profile 
obtained numerically from Eq. 10 compares very well 
with the actual computer simulation results of the ran- 
dom lattice model (see Fig. 01, as expected. 
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IV. DISCUSSION 

In this paper, we have extensively studied structure 
formation and their dynamics in a one-dimensional in- 
elastic lattice model for granular gases. These structures 
appear in the form of clusters separated by shocks. The 
locations of the shocks in this lattice model at early times 
are decided from the initial configurations of the parti- 
cles' velocities in a very robust manner. In order to pre- 
dict the locations of the shocks, we have demonstrated 
a procedure to process the initial velocity configurations 
data of the particles, and this procedure works eq ually 
well also for the one-dimensional inelastic gas [lj, fl7| . 
The coefficient of restitution does play a role to decide 
which of the large jumps in the particles' velocities at 
early times yield shocks, but at late times, the macro- 
scopic velocity profile for a given model with a given 
initial configuration of the particles' velocities is inde- 



pendent of the coefficient of restitution. 

Thus, we observe that in terms of detailed structure 
formation, systems with any r < 1 "flow" towards the 
sticky limit r — 0. Such "universality" has been found in 
Ref. |17| in terms of global quantities, while our results 
suggest that a broader "universality" holds for micro- 
scopic quantities, such as the locations of shocks (when 
the same initial velocity configuration is considered). 

In addition, for the random lattice model, we have also 
studied the dynamics of an isolated cluster and cluster- 
cluster collisions in detail and obtained an effective "hy- 
drodynamic" equation. We hope that the analyses pre- 
sented here can be successfully used to study coarsening 
problems in realistic granular gases (i.e., in two or three 
dimensions). 

S.O. and D.P. are financially supported by the Dutch 
research organization FOM (Fundamenteel Onderzoek 
der Materie). 
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